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Abstract 

We present the results of a search for neutral Higgs bosons produced in a- 
ssociation with b quarks in pp — > bbip — » bbbb final states with 91 ± 7 pb _1 of 
pp collisions at = 1.8 TeV recorded by the Collider Detector at Fermilab. 
We find no evidence of such a signal and the data is interpreted in the context 
of the neutral Higgs sector of the Minimal Supersymmetric extension of the 
Standard Model. With basic parameter choices for the supersymmetric scale 
and the stop quark mixing, we derive 95% C.L. lower mass limits for neutral 
Higgs bosons for tan j3 values in excess of 35. 
PACS numbers: 13.85.Rm, 12.15.Ji, 14.80.Cp, 14.80.Bn 
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A fundamental question which remains open today in particle physics is the origin of the 
electroweak symmetry breaking. The simplest mechanism in the Standard Model (SM) and 
in many supersymmetric extensions is spontaneous symmetry breaking, achieved through 
the introduction of one or more scalar field doublets. The SM assumes one doublet of scalar 
fields and a single physical Higgs boson (hsu), with unknown mass but with fixed couplings 
to other particles. A more complex symmetry breaking mechanism occurs in the Minimal 
Supersymmetric extension of the Standard Model (MSSM), where several physical scalar 
states are predicted: three neutral bosons (the CP-even h and H, and the CP-odd A) and 
two charged bosons (H^). A distinct feature of the MSSM is the modified couplings of the 
Higgs particles, in particular the enhancement of the bottom-Higgs Yukawa couplings by 
tan p for the case of the bbA vertex, where tan /3 is the ratio of the vacuum expectation 
values of the two Higgs doublets of the theory. The Higgs sector of the MSSM is completely 
determined at tree level by two free parameters, chosen to be the mass of the CP-odd Higgs 
boson, rriA, and tan/?. The mass of the CP-even Higgs boson h, rrih, is constrained to be 
less than mzo \ cos2/3|. Radiative corrections substantially modify the masses and couplings 
of the two neutral CP-even Higgs scalars, in particular the upper bound on m h [|TJ. 

At the Tevatron, one of the Higgs production mechanisms likely to be observed in both 
the SM and some regions of the MSSM parameter space is the associated production pp — > 
V + ip, where V = W, Z and <p = h,H, hsM- CDF has already reported on searches for this 
channel (V + Ksm) with different signatures |||3|]. In this Letter we exploit the enhanced 
bottom-Higgs Yukawa couplings of the MSSM to test the large tan j3 sector of the theory by 
searching for the process pp — > bbip — > bbbb with ip = h,H, A. Our sensitivity in this search 
is limited to the region of parameter space corresponding to tan/? > 35. In this region, at 
least one of the CP-even Higgs bosons has similar couplings to, and is always degenerate 
with, the CP-odd Higgs boson. For itla < 110 — 125 GeV/c 2 (depending on tan /3 and on the 
parameters of the stop quark mass matrix), mh — ttia, while for rriA ~ HO — 125 GeV/c 2 , 
rriH — ttia- Therefore this analysis covers a simultaneous search for two or more Higgs 
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signals with an experimental signature of four b jets in the final state. 

The search reported here is based on 91 ±7 pb _1 of integrated luminosity recorded during 
the 1994-95 Tevatron run. The CDF detector is described in detail elsewhere The 
silicon vertex detector (SVX) consists of four layers of axial microstrips located immediately 
outside the beampipe with an innermost radius of 2.9 cm ||. The SVX provides precise 
track reconstruction in the plane perpendicular to the beam and the ability to identify 
secondary vertices produced by heavy flavor decays. The momenta of charged particles are 
measured in the central tracking chamber (CTC), which lies inside a 1.4 T superconducting 
solenoidal magnet. Outside the CTC, electromagnetic and hadronic calorimeters arranged 
in a projective tower geometry cover the pseudorapidity region \r)\ < 4.2 and are used 
to identify jets. The data sample was recorded with a trigger which requires four or more 
clusters of contiguous calorimeter towers, each with transverse energy Ex > 15 GeV, and a 
total transverse energy J2 Et > 125 GeV. 

The initial steps of the data selection are the same as in the recent CDF SM Higgs 
search ||. We start by rejecting cosmic ray events, beam halo, and detector noise. Events 
with one or more identified electrons or muons from vector boson decays defined as in |2j 
are also rejected. After this selection, events are required to have at least four jets with 
Et > 15 GeV and well contained within the fiducial calorimeter regions \q\ < 1.5. Jets 
are defined as localized energy depositions in the calorimeters and are reconstructed using 
an iterative clustering algorithm with a fixed cone of radius AR = y/ Arj 2 + A<p 2 = 0.4 in 
r] — <p space [0. Jet energies are corrected for energy losses in uninstrumented detector 
regions, energy falling outside the clustering cone, contributions from underlying event and 
multiple interactions, and calorimeter nonlinearities. The selected data sample is dominated 
by QCD multijet events and contains 207, 604 events. The four highest--Er jets in an event 
are then ordered in Et and a search sample is obtained for each of a set of Higgs boson 
masses by requiring the three highest--Er jets to pass cuts which are Higgs boson mass 
dependent. This is motivated by the fact that the Et spectrum of the leading jets for the 
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signal is, on average, larger than the QCD background, and grows with increasing scalar 
boson mass (m v ). We use the leading order (LO) parton level matrix elements || encoded 
in the PYTHIA v5.6 Monte Carlo program || along with a full simulation of the CDF 
detector to simulate the signal. We use the CTEQ3L parton distribution functions and 
set a factorization scale equal to the Higgs boson mass in the simulation. We find optimal 
Et cuts by maximizing the expected significance of the signal. For a Higgs boson mass of 
m^p = 120 GeV/c 2 these cuts correspond to 48, 34 and 15 GeV for the leading jet, second 
and third leading jets of the event, respectively, and vary roughly linearly with the Higgs 
boson mass. We then require that at least three among the four highest- Et jets in the event 
are identified (tagged) as b quark candidates. The algorithm used to identify secondary 



vertices |10| begins by searching for jets which contain three or more displaced tracks. If 
none are found, the algorithm searches for two-track vertices using more stringent track 
criteria. A jet is tagged if the transverse displacement of the secondary vertex from the 
primary vertex exceeds three times its uncertainty. A requirement on the azimuthal angular 
distribution of the two highest- E T 6-tagged jets in the event, Aip b i > 109°, reduces the heavy 
flavor QCD content of the sample attributed to gluon splitting. This cut preserves ~ 90% of 
the signal events which favor a larger angular separation between the 6-tagged jets coming 
from the Higgs decay. After the three tag requirement and the A(p cut we are left with 20 
and 13 events, respectively, for the case of the m v = 70 GeV/c 2 selection. 

To reconstruct the signal we select one of the possible invariant mass combinations bet- 
ween the highest-.&r jets in the event. From Monte Carlo we find that the mass of the 
highest-Er jets in the event (m 12 ) for signal masses m v > 120 GeV/c 2 , and the mass of the 
second and third highest--&r jets in the event (m 2 3) for signal masses < 120 GeV/c 2 
enhance the signal resolution, ^(m^/m^. The use of these distributions also minimizes the 
percentage of events for which at least one of the jets in the dijet mass is not associated 
with a b quark from a Higgs boson decay. All signal mass distributions contain a Gaussian 
core with a resolution which depends on the Higgs boson mass and varies from ~25% for 
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rritp < 120 GeV/c 2 to ~13% for m v > 120 GeV/c 2 . We increase the expected significance 
of the signal by applying mass window cuts on the m 12 and m 23 distributions which vary 
between ±15(m^) and ±35 (m v ), depending on the Higgs mass, and centered on the mean 
of the fit distributions. A cut on the invariant mass distribution between the two 6-tagged 
highest- .&r jets of the event further discriminates against heavy flavor QCD events. All 
mass cuts were chosen to maximize the expected significance of the signal. Table | shows 
the number of observed triple 6-tagged events left after all cuts as a function of mass. Five 
events are left after all cuts in the mass bin at 70 GeV/c 2 . All events for the mass bins 
above 70 GeV/c 2 are included in this sample of five events. 

In addition to the large QCD multijet background, other sources of heavy flavor in the 
triple 6-tagged sample include multijet tt production (t — > Wb, W — > qq'), Wbb and Wcc 
with W — > qq', Zbb and Zee with Z — > bb/cc, and fake triple-tags. They are estimated 
from a combination of Monte Carlo and data. We use the HERWIG v5.6 Monte Carlo 
generator JO] with the CDF measured cross section [a t i = 6.5^4 pb) |12| and a top mass 
of m t = 175 GeV/c 2 to predict the expected number of tt events. Electroweak processes 
are also estimated with the same Monte Carlo generator. Fake triple-tags are defined as 
events in which at least one of the three 6-tagged jets contains a false secondary vertex in 
a light quark or gluon jet. Fake tag probabilities are parameterized by measuring in several 
inclusive jet data samples the fraction of jets in which a secondary vertex is reconstructed 



on the wrong side of the primary vertex with respect to the jet direction [10.13 



The same fit technique that was used to estimate the QCD heavy flavor normalization 
in the SM Higgs search || is applied to this analysis. We reduce the 6-tag cuts on our data 
sample from triple to double 6-tags. This gives a high statistics background-rich sample 
in which we fit the double 6-tagged dijet mass distribution to a combination of signal and 
SM backgrounds. The shape of the QCD heavy flavor distribution is obtained from the 
PYTHIA Monte Carlo program. We generate all QCD jet production channels, and retain 
the events that contain a heavy quark produced either in the hard scattering or in the 
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associated radiation process. The signal and QCD normalization is left free in the fit while 
the SM non-QCD background, both shape and normalization, is obtained from Monte Carlo 
simulation. The QCD normalization in our triple 6-tagged sample is then obtained from 
the ratio of normalizations of double 6-tagged to triple 6-tagged QCD events taken from the 
PYTHIA Monte Carlo simulation. Table |TJ lists the individual QCD, ti, Wbb, Wee, Zbb, 
Zee, fake triple tags, and total expected contributions to the final observed sample as a 
function of mass. From these numbers we find no evidence for the presence of a Higgs boson 
signal. Figure [l] shows the myi and 77223 distributions for the observed triple 6-tagged sample 
compared to the SM background expectations and for three different selections corresponding 
to rriip = 70, 120, and 200 GeV/c 2 . Also shown are the signal plus background shapes 
normalized to the expected number of events for tan/3 = 50 and the case of no mixing in 
the scalar top sector (no mixing scenario). 

We calculate the signal detection efficiencies and normalizations from Monte Carlo. The 
overall detection acceptances with their total uncertainties are shown in Table | as a function 
of mass. They are within a range of 0.2% to 0.6%, increasing with mass. This low acceptance 
is dominated by the small multijet trigger efficiencies (~1% to ~7%, increasing with the 
signal mass) and, to a lesser extent, by the triple 6-tag requirement (~20%). The low values 
for the trigger efficiency are due to the high Et thresholds and multiplicity requirements 
on jets. The trigger efficiency curves have been obtained with a trigger simulation with 
parametrized curves estimated from data. The total systematic error includes uncertainties 
in the multijet trigger simulation (5%), in the modelling of gluon radiation (10% to 7%, 
depending on the mass), in the calorimeter energy scale (10% to 2%, depending on the 
mass), in the luminosity measurement (7%), and in the 6-tag efficiencies (10%). 

From the data in Table | we set upper limits on bbip — > bbbb (if = h, H, A) production 
using a one-sided conditional frequentist construction ]TJJ, where systematic uncertainties 
are approximately taken into account by Bayesian averaging over the systematic parame- 
ters assuming gaussian a-priori distributions around their best estimates. The 95% 
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CL. upper limits on the total expected signal events as well as on the production cross 
section times branching fraction are listed in Table |. Using the LO theoretical cross sec- 
tions for o~(pp — ► bb(p)BlZ((p — > bb) with (p = h, H, A and the bottom-Higgs Yukawa coupling 
calculated with a running bottom quark mass evaluated at the Higgs boson mass scale, 
TUbirn^) ~ 3 GeV/ c 2 , we exclude regions of parameter space for m h — tan (3 and — tan j3, 
respectively, as shown in Figures ^| and [| Results are shown for two common choices 
of the stop quark mixing parameter Q: no mixing (A t = /icot/3), and maximal mixing 
(A t = \i cot f3 + VQms, with \x the supersymmetric Higgs boson mass parameter and A t a 
soft breaking parameter). In all cases we set ms, the supersymmetric mass scale, to be 1 
TeV/c 2 and m t = 175 GeV/c 2 . As a test of the sensitivity of our calculated limits to the 
background estimate, if we assume a zero background hypothesis the result limits increase, 
for all signal masses, by 5(tan/3) = 10. 

In conclusion, we have searched for neutral Higgs bosons produced in association with b 
quarks through the reaction pp — ► bbip — ► bbbb. We do not find evidence for the presence of 
a signal and 95% CL. upper limits are set on the production cross section times branching 
ratio. The results have been interpreted in the context of the MSSM Higgs sector and lower 
mass limits for neutral Higgs bosons derived for tan/3 values in excess of 35. 
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TABLES 



TABLE I. Number of observed and expected background events after the final selection and 
for the different SM contributions as a function of mass. Last three columns show the total 
acceptances and expected 95% C.L. upper limits on the number of signal events (N signa i) and on 
a x BIZ, respectively. 



rriip 


Observed 


Expected 


Acceptance N signa i 


a x BIZ 


(GeV/c 2 ) 


Events 


Background 


(%) 


(95% C.L.) (pb, 95% C.L.) 


7D 


o 


4 fi + 1 4 

"4:.U H_ J. . t: 


0.16 ±0.03 7.9 




53.3 


80 


4 


4.6 ±1.4 


0.22 ±0.04 6.6 




31.7 


90 


3 


3.8 ±1.3 


0.23 ±0.04 5.8 




27.1 


100 


3 


3.8 ±1.3 


0.25 ±0.04 5.9 




25.7 


110 


2 


3.7 ± 1.1 


0.25 ±0.04 4.8 




20.7 


120 


2 


3.5 ±1.1 


0.28 ±0.05 4.9 




19.2 


130 


1 


2.6 ±0.9 


0.28 ±0.05 4.1 




15.8 


140 


1 


1.7 ±0.8 


0.29 ±0.05 4.3 




16.2 


150 





1.5 ±0.8 


0.30 ±0.05 3.2 




11.5 


200 





1.2 ±0.7 


0.41 ±0.07 3.2 




8.5 


250 





1.0 ±0.7 


0.47 ±0.08 3.2 




7.5 


300 





0.1 ±0.4 


0.59 ±0.09 3.2 




5.8 


TABLE II. Expected QCD, fake triple tags, tt, Wbb + Wcc, Zbb + Zee, and total number of 


background events as a function of mass. 










m v (GeV/c 2 ) 


QCD 


Fakes 


tt 


Wbb + Wcc 


Zbb + Zee 


Total 


70 


2.97 ±0.70 


0.5 ±1.2 


0.70 ±0.18 


0.09 ± 0.06 


0.37 ±0.02 


4.6 ± 1.4 


80 


2.97 ±0.70 


0.5 ±1.2 


0.70 ±0.18 


0.09 ± 0.06 


0.37 ±0.02 


4.6 ± 1.4 


90 


2.16 ±0.55 


0.5 ±1.2 


0.70 ±0.18 


0.09 ± 0.06 


0.37 ±0.02 


3.8 ± 1.3 


100 


2.16 ±0.55 


0.5 ±1.2 


0.70 ±0.18 


0.09 ± 0.07 


0.37 ±0.02 


3.8 ± 1.3 


110 


2.16 ±0.55 


0.4 ±0.9 


0.68 ±0.18 


0.09 ± 0.07 


0.37 ±0.02 


3.7 ± 1.1 


120 


2.16 ±0.55 


0.3 ±0.9 


0.66 ±0.17 


0.07 ±0.06 


0.29 ±0.02 


3.5 ± 1.1 


130 


1.44 ±0.45 


0.3 ±0.8 


0.64 ±0.17 


0.05 ±0.05 


0.21 ±0.02 


2.6 ±0.9 


140 


0.73 ± 0.40 


0.2 ±0.7 


0.60 ±0.16 


0.04 ±0.05 


0.16 ±0.02 


1.7 ±0.8 


150 


0.73 ± 0.40 


0.1 ±0.7 


0.55 ±0.15 


0.03 ± 0.05 


0.12 ±0.02 


1.5 ±0.8 


200 


0.73 ± 0.40 


0.1 ±0.6 


0.32 ± 0.09 


0.00 ± 0.00 


0.00 ±0.00 


1.2 ±0.7 


250 


0.73 ± 0.40 


0.1 ±0.6 


0.17 ±0.05 


0.00 ± 0.00 


0.00 ±0.00 


1.0 ±0.7 


300 


0.01 ±0.01 


0.1 ±0.4 


0.00 ± 0.00 


0.00 ± 0.00 


0.00 ±0.00 


0.1 ±0.4 
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FIGURES 
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FIG. 1. Invariant mass distributions m\2 and TO23 for the observed triple 6-tagged sample co- 
rresponding to the = 70 GeV/c 2 (top), = 120 GeV/c 2 (middle), and = 200 GeV/c 2 
(bottom) selections. The data is compared to the expected QCD only background, the total SM 
backgrounds, and the total background plus signal for tan/3 = 50 and the no mixing case. The use 
of for < 120 GeV/c 2 , and m\2 for > 120 GeV/c 2 maximizes the fraction of correct jet 
assignments and enhance the signal resolution (see text). The mass cuts are not applied. 
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(a) m,(GeV/c 2 ) ( b ) m h (GeV/c 2 ) 



FIG. 2. CDF 95% C.L. excluded region in the parameter space m^ — tan/3 for the two stop mix- 
ing scenarios: (a) no mixing, and (b) maximal mixing. Also shown are the theoretically forbidden 
regions and the LEP exclusion region for their no mixing and rn™ ax scenarios [16]. 
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FIG. 3. CDF 95% C.L. excluded region in the parameter space — tan/3 for the two stop 
mixing scenarios: no mixing (dashed lines) and maximal mixing (solid line). Also shown is the 
LEP exclusion region for the no mixing scenario [16]. 
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